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ABSTRACT 


A design model for a once-through waste heat recovery 
unit with a segmented fin-tube arrangement was developed 
along with a simple model of a combined gas and steam 
(COGAS) turbine propulsion system. These models were inte- 
grated and applied in a computer program written in FORTRAN 
IV for the IBM 360-67 computer. Waste heat recovery unit 
designs were produced and tested at off-design conditions. 
Using the space constraints and power requirements of a 
Navy destroyer-type ship, one design was selected and 
employed to make estimates of possible fuel savings to be 


realized through the application of a COGAS system. 
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i. INTRODUCTION 


A. BACKGROUND 

The U. S. Navy has made a commitment to gas turbine 
propulsion systems, with two classes of gas turbine ships 
entering the fleet (DD-963 and FFG-7) and one class pro- 
posed (DDG-47) at the time of this writing. The DD-963 class 
has two shafts with two General Electric LM 2500 gas tur- 
bines powering each shaft. The FFG-7 class is a single- 
Shaft ship with two LM 2500 gas turbines. The DDG-47 class 
propulsion plant, as proposed, will be similar to that of 
the DD-963 class. 

Currently, means for conserving fuel are being sought 
in all sectors of the economy and government. In most of 
our steam turbine-powered ships fuel conservation must be 
pursued largely through operating practices such as reduced 
steaming time or reduced speed when steaming. Several 
alternatives exist, however, for fuel consumption reduction 
for the gas turbine propulsion system. These alternatives 
involve combining another type of propulsion system with 
tne gas turbine in order to reduce the overall fuel consump- 
tion by either operating the combined systems in parallel 
or by operating the less expensive system in the cruise 
mode. Some examples of these combinations are: 

Gas Turbine and Diesel (CODOG, GODAG) 
Gas Turbine and Steam (COGAS) (steam system not 


separately fired) 
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Gas Turbine and Steam (COSAG) (steam system 
separately fired) 
The COmbined Gas turbine And Steam turbine (COGAS) system 
is probably the best of these alternatives with respect to 
Buel CONSUMpPtion, initial cost, maintenance, and manning 
requirements. 

The COGAS ship propulsion plant uses the exhaust heat 
from a gas turbine to generate steam in a waste heat re- 
covery unit (WHRU). This steam is used to drive a steam 
turbine which is operated in parallel with the gas turbine, 
for ship propulsion. A common reduction gear connects both 
turbines to the propeller shaft. Figure 1 illustrates the 
COGAS system. 

Several studies have examined the potential of the COGAS 
system for improving propulsion plant fuel efficiency. 
Giblon and Rolih [Ref. 1] studied the commercial ship 
application of COGAS. References 2 and 3 investigated the 
feasibility of the COGAS system for military ship application. 
In all three studies the COGAS system was designed around 
a gas turbine brake horsepower of about 23,000. Additionally, 
all three employed a drum-type boiler as the waste heat 
recovery unit. 

In this study, COGAS system designs were considered for 
a destroyer-sized ship with two shafts and two LM-2500 gas 
turbines for each shaft. One gas turbine was equipped with 
a WHRU. In the COGAS mode of operation the WHRU-equipped 


gas turbine was operated in parallel with the steam turbine, 
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powering one shaft. The other shaft was allowed to drag. 

It was assumed that the principal employment of the COGAS 
system would be for "cruise" conditions. That is, the COGAS 
system would be operated primarily for steady steaming. The 
horsepower range over which the COGAS system could be 
operated is, theoretically, limited only to the range of 
powers over which the gas turbine can be operated. Thus, 
the power range would vary from the COGAS system power 
Output with the gas turbine at idle to the system power 
output attained with the gas turbine operating at its maxi- 
mum continuous power rating. 

When the COGAS system operating range is selected, the 
desired goal of fuel conservation must be balanced against 
system reliability and routine maintenance requirements. 

If only one gas turbine is equipped for COGAS operations, 
the system downtime will increase as the turbine is operated 
at higher power settings. In order to increase the time 
between major gas turbine maintenance requirements, two 
COGAS systems could be installed in a twin-shaft ship. The 
increased overall reliability which this redundancy would 
afford, however, would be at least partially offset by the 
higher initial cost, the increased manning requirements, 

and the added weight of two systems. 

A once-through counter-cross-flow heat exchanger was 
Selected as the WHRU. Reference 6 states some of the 
principal advantages of the once-through steam generator 


Over the drum-type boiler. These advantages are as follows: 
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lower initial cost, faster response, more compactness, 
and lower operating cost. Additiorfally, the once-through 
unit provides a somewhat simpler model than the drum-type 
heat exchanger. A schematic representation of a once- 
through WHRU installed in the exhaust ducting of a gas 


turbine is shown in figure 2. 


fe «OBJECTIVES 

There were five main objectives for this thesis. 
1. Develop a design model for the once-through WHRU. This 
model waS to size a WHRU and predict is performance given 
initial conditions such as feedwater inlet temperature, 
steam outlet temperature, WHRU operating pressure, inlet 
and outlet gas conditions, and gas flow rate. 
2. Consider several sets of initial conditions and design 
constraints (e.g., maximum frontal dimensions) and produce 
WHRU designs for each set. 
3. Evaluate these designs and integrate them with a simple 
COGAS system model in order to test the design model and 
to develop a clearer understanding of the WHRU performance 
in the COGAS system. 
4. From the evaluated designs, predict the performance 
of the COGAS system. 
5. Provide a framework for future studies involving a 
detailed, complete COGAS system cycle analysis and 


Optimization. 
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IT. MODEL DESCRIPTION 


A. OVERVIEW AND INITIAL CALCULATIONS 

The waste heat recovery unit (WHRU) is modeled as a 
once-through counter-cross-flow heat exchanger. The model 
is applied in a computer program for use in heat exchanger 
design for a combined gas and steam (COGAS) turbine pro- 
pulsion plant system. The same model could be adapted 
for use in a complete cycle analysis of the COGAS system. 

The WHRU model is divided into three principal sections: 
heating, boiling and superheating. The specified initial 
conditions for the model are: (1) gas temperature into 
the heat exchanger (Toy) (2) minimum average gas exit 
temperature from the heat exchanger (Tyg), (3) gas flow 
rate (mo), (4) water inlet temperature (Te), (5) water/ 
steam pressure (P-). Useful background material for the 
model formulation was obtained from references 1-5. 

Two simplifications should be noted at this point. First, 
the water-side pressure drop in the WHRU is neglected. 
Second, gas temperatures, evaluated at several points in 
the model, are average temperatures. Specifically, the 
average inlet and outlet gas temperatures for the heating 
and superheating sections and the average inlet and outlet 
gas temperatures for each pass in the boiling section are 
evaluated. In each of these calculations the gas tempera- 


ture distribution is neglected. Clearly, even if the gas 
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turbine exhaust flow is uniform in temperature as it 
approaches the WHRU, the gas temperature will not be uni- 
form after the first WHRU pass is encountered. Due to the 
varying gas-to-fluid temperature difference across each 
pass the gas temperature will be distributed non-uniformly 
after the first WHRU pass. One can, however, calculate the 
average gas temperature approaching a particular pass 

based on the heat release from the gas in the previous pass 
and the average gas temperature entering the previous pass. 
Averaged gas temperatures calculated in this manner are 
used throughout the model. 

Prior to the heating section, some initial calculations 
are performed to determine the water/steam mass flow rate 
and to set the pinch point. In order to calculate the 
water/steam mass flow rate an energy balance is performed, 
uSing the starting conditions, as follows: 

Q = ma, Se Be Sede. 
where C g 1s evaluated at the average gas temperature 


= (T Ne) ee The fluid mass flow rate 1S 


Toavg eal 


calculated from 


where Ney and Neg are inlet and outlet fluid enthalpy 
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respectively. Intermediate fluid and gas temperatures are 
calculated as follows: 
1. The heating section is assumed to heat the water from 


inlet conditions to the saturation temperature, A 


Te5° 
heat balance is performed on the water side yielding a 


heating section heat transfer of 


The required average gas temperature entering the heating 


section is then calculated from 


2. In the boiling section it 1S assumed that sufficient 
heat 1s added to produce saturated vapor at the outlet. 
An energy balance on the steam-Side provides the total heat 


transfer in the boiling section, 


where Ae 1s the enthalpy of saturated vapor at the 
pressure specified. The required average gas temperature 


entering the boiling section is 


Ze 





3. In the superheating section, the steam temperature is 
raised from the saturation temperature to the outlet 
temperature specified in the initial conditions. The heat 


transfer in the superheating sections is 


where Neg is the enthalpy of steam corresponding to the 


steam outlet temperature T The gas inlet temperature 


£4° 
aa is known from the initial conditions. 

imew pinch point”, for the purposes of this model, is 
defined as the difference between the gas temperature enter- 
ing the heating section and the water temperature leaving the 
heating section (saturation temperature). This temperature 
difference will normally correspond to the smallest tempera- 
ture difference between gas and fluid in the heat exchanger. 
The only other likely location for the smallest temperature 
difference between gas and fluid is at the gas inlet, or 
Superheater outlet. Figure 3 depicts the typical gas and 
fluid temperature distribution in the WHRU. 

The model must provide some internal control over the 
Pinch point temperature difference as defined above. This 
1s desirable because the gas-fluid temperature difference 
at the fluid saturation point is indicative of the amount of 
heat transfer area which will be necessary to heat the in- 


coming water to the saturation point. As the temperature 


26 





TYPICAL WHRU GAS-FLUID TEMPERATURE DISTRIBUTION 
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difference at the pinch point increases the heating portion 
of the heat exchanger becomes more efficient from the point 
of view of heat transfer area required. The minimum pinch 
point temperature difference is set at 25 F. After the 
interim temperatures are established the pinch point temp- 
erature difference 1s checked. If the difference is less 
than 25 F, the gas temperature at the heating section inlet 


is reset to 


From the overall energy balance 


( a) — Tl 


(eS Teas 
exe) ey ene g4 


and the energy balance across the heating section 


Ca ec TS = m_(h_-.-h : 
pg g | g3 g4! 


a revised gas outlet temperature may now be calculated from 


= 3 il 
54 l-a 
where 
_ ee Bae, 
es pe er 
£4 sm 
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After this new gas inlet temperature is established, all 
initial calculations are performed again, to establish a 
new water/steam mass flow rate and to recheck the pinch 
Bount temperature difference. Once these calculations are 
complete, the geometric parameters for the heat exchanger 


must be specified. 


eee GEOMETRY 

A fin-tube with helically curved extended surfaces on 
circular tubes is assumed for this model. The fins are 
segmented. The tubes are configured in banks of one row 
each and the rows are staggered. The base tube surface 1s 
taken from Ref. 7. The description of the finned tubes 


is as follows. 


dq; = tube inside diameter = 1.86 in. 

dq = tube outside diameter = 2.00 in. 

dq = fin root diameter = 2.00 in. 

Ne = fins per inch = 5.94 

QR Sec nervagnt = 1.015 in. 

Re SP onGiiror ut trom flinetip =. .0.82 in. 

d = fin outside diameter = 4.03 in. 

a Ss fin thickness = 0.048 in. 

Wew cs fin segment width = 0.17 in. 

er on number of segments in 360 degrees = 38 
The finned tube configuration is shown in figure 4. The 


P@eemtength, L, and number of tubes per row Me are 


chosen by the designer according to the frontal area, A,, 


Zo 





Ws 


f~ 


s 
“nee | 


Figure 4: Segmented Fin-Tube Configuration 


desired. The tube layout is shown below in figure 5. 

The center-to-center tube spacing in the transverse direc- 
tion is 4.50 inches. The equilateral, staggered tube 
arrangement used in this model leads to a spacing normal 
to the gas flow, Sal of 4.50 inches and a spacing parallel 
to the gas flow, per of 3.90 inches as shown in figure 5. 
The heat exchanger height is established by the number 

Of WHRU passes (rows) and the tube spacing parallel to the 


gas flow se) ¢ 
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Figure 5: Model Tube Layout 


In order to establish the minimum gas flow cross- 
sectional area the total "blocked" frontal area, Abr OL 


the heat exchanger must be calculated from 
is 


The total inside area available for heat transfer per pass 


1s 
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The gas side surface area available per pass for heat 
transfer is the sum of the fin surface area and the bare 
tube Surface area per pass. The fin surface area per tube 


is calculated from 


Acin = Ne LIN. (2 de wo +2 te, eae Wee te) 
Nee 2 2 
= x (dep ae 
where dey = d. = 2k. The bare tube area per tube is 
Any = 7 ae Ly pp dq te Ne L 


Therefore, the total outside area per pass available for 


heat transfer is 


Boo = Nesp ‘Atin oe 


The cross-sectional area for fluid flow is calculated from 


EP ee 
Bee = 7 a Ae 
With the mass flow rates, terminal temperatures and heat 
exchanger geometry established, the remainder of the model 
may be solved for number of passes, actual interim tempera- 


tures, location of the fluid phase changes, and the gas 


Side pressure drop. 


Sy 





C. HEATING SECTION 

The gas~side Reynolds number for the heating section is 
calculated initially using the gas bulk temperature to find 
the gas properties. With this Reynolds number, a j-factor 
Mme friction factor, £, are obtained from polynomials fit 
to the data for tube layout number 5 in Ref. 7 (fig. 5). 
The friction factor will be used later to calculate the 
gas side pressure drop in the heating section. The j-factor 
is related to the heat transfer coefficient Ae by the 


following relationship. 


By introducing 





N N 
{= u p2/3 ie on : 
R ne y : 
eg xg Rag 2 rg 
and 
Neeser Pal 
u 1S | 
where 
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Therefore, a relationship may be written for the heat 
transfer coefficient as follows: 
K 
—— ze R ee 
eg rg 
O 
The water-side heat transfer coefficient is calculated by 


using the Dittus-Boelter correlation 


where all properties are obtained at the bulk temperature 

of the water in the heating section. Using the tube wall 
resistance together with Ne and Nes the overall heat transfer 
coefficient for the heating section may be written in terms 


of the inside area as 





iG , a 
‘ape l i Ba in(d_ ) ae l 
a. 27K, Ney L La ny h, 
where 
= ” _ Poet 
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The fin efficiency, Nee 1s calculated from the expression 


tanh ML 
Ne Mee 


where 


a Vh P/KA 


and, if the fin 1S approximated by a set of rectangular 


mo 


2 


the cross-sectional area, A, of a rectangular strip may 


Strips extending from the tube wall, L = Now, 


be written as 


So, uSing known geometric parameters and the thermal conduc- 
tivity of the fin metal, the parameter ML may be expressed 
as 


ML = C vh_ 
g 


where 
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Now, 


tanh C vh Acsn 
n. = (L- (2 - ———— 4) >—= 
€ oF Op 


The number of passes required inthe heating section is 
calculated, using the effectiveness-NTU method, in the 
following way. 

1. An average pass effectiveness for the heating section is 
calculated from the following expression for cross-flow 
effectiveness with both fluids unmixed obtained from 


Ref. 8, 


exp(“NCn) - 1) 


where 


p Cmin pé Me 
See PG G 
ves 
N = NTU = == : 
Min 
: eaoe22 


Water and gas properties are taken at bulk temperatures. 
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2. An overall heating Section effectiveness is calculated 


from 


i, speek ere, oC 
p min max) Nn ea 
f=. 
oa =" ele Cre 7% Coe: 
Pp min max )n _ min 
le-ece (S 
p max 


where n = number of passes. This expresSion was obtained 
from Ref. 9. The overall effectiveness formulation was 
actually derived with the condition that fluids are mixed 
between passes. It can be shown [Ref. 9], however, that 
the error is not large when the expression is used for the 
case where fluids are unmixed between passes. 

3. The expression for ae 1s now used in an iterative 
fashion starting with n = number of passes = 1. Each time 

E Oa 1s solved, the gas temperature into the heating section, 


fee tS found from 
aS 


E Caras, L aS fi 
s: oa min gl max g4 
Ab = 
g3 E Cee 
ee) Hue se max 


Poreh is derived from 


Smase 29 7 Tea? 


E ” oa ae 
Oa Catia ae ea 


4. Once the inlet gas temperature is established for a 


particular number of passes, an energy balance may be 


Sah 


performed to solve for the total heat transfer which would 
take place in the heating section for that number of passes. 


rane dls ae 
h g “pg '"g3 g4? 
Ber fhe enthalpy of the water at the outlet of the heating 


section for a particular iteration may be found from 


and this yields the temperature of the water at the outlet, 


a for the number of passes under consideration. 


inl 
6. The outlet water temperature, Teo: may be checked to 

see if it exceeds the saturation temperature. 

7. This set of calculations is performed until Te5 exceeds 
the saturation temperature. This naturally means that 
boiling begins in the last pass of the current total number 
of passes. Initiation of boiling wlll be treated in the 
boiling section of the model. Therefore, the heating section 
ends with the pass prior to that which initiates boiling. 

The results of the above calculations are: (a) The total 
number of passes contained in the heating section, (b) the 
new gas temperature at the inlet to the heating section Beat 
and (c) the new water outlet temperature Te: 
With these new heating section terminal temperatures an 


average outside wall temperature for the section may be 


calculated from 
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where A, ;/AL are total inside/outside areas for the section 
and Tob/ eb are the new gas/fluid bulk temperatures for the 
section. The expression for eee is derived from the follow- 


ing formulation for heat transfer in the heating section, 





_ Tgp 7 Tp _ gb 7 Two 
ae R 
YR O 
which reduces to 
Ss 
wo | Tob 7 YR ( gb 7 Tey) 
2 ical! 
where 
L 
R —_ 
Oo Ny ey Ao 
and 
= JL, 
Rep aa 
on tL 
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This gas-side film temperature 1s now introduced at the 
beginning of the calculations for the heating section, to 
replace the gas bulk temperature. All calculations are 
performed again with Tye: 

The gas-Side pressure drop in the heating section may 
now be calculated using the previously obtained friction 
factor, f, from the correlation contained in Ref. 7. The 


gas-side pressure drop is 


2 


AP _ NE See Ne/t Mw) 0.14 
g p Uy 
where 

m 
G = SS 
max A_. 

min 
oe = number of tubes per row 
fe) = gas density 
Le = gas viscosity at wall temperature 
Wy = gas viscosity at bulk temperature. 


Eee SOLILING SECTION 

The boiling section is solved pass-by-pass since most 
of the section will involve two-phase flow on the cold side 
and the inside heat transfer coefficient will change with 
amettty, x, and heat flux, q". The correlation for the 


two-phase region inside heat transfer coefficient selected 
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for this model is that recommended by Tong [Ref. 10] for 
both nucleate boiling and forced convection. The equation 
for the two-phase flow heat transfer coefficient was 


given by Schrock and Grossman in Tong [Ref. 10] as 





h 
Tt n 
A - Slew + ag). 
Q fg Pe 
The constants are given by Wright [Ref. 10]: B= STOO 


A= me 10, n= 0.66. The heat transfer coefficient 


assuming a total liquid flow is 


eet a) 08 oe “2.0.1 
got +} — 198 BE 41°? 
i 2 2 


C 


Ay = 0.023 


and the Martinelli parameter is 





Maen limits Of the data for this correlation are 


Milaieleys we. 05 -— 0.57 


2) 
IMEC Ie IAN eH 350 20" = Ue crenttoye edge 5 ey 


In this model, heat flux in the boiling section will be at 
levels considerably below eos BTU/hr-£t* and the full 
range of quality must be considered. Therefore, in order 
to observe the performance of the Schrock and Grossman 
correlation at lower heat flux and quality above 0.57 a 


plot was made (fig. 6) of the ratio he pe/ Bp vs. quality 
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Dengler & Addoms 


Ne pe Ay vs. QUALITY 


PEGE vO 
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where hy 1s the heat transfer coefficient from the Dittus- 
Boelter correlation for water at saturated liquid condi- 
tions. The correlation of Dengler and Addoms as given in 
Tong [Ref. 10] for entirely forced convection vaporization 
mormecluded for comparison. Tong [Ref. 10] indicates that 
low heat flux and high water-vapor mixture velocity favor 
the forced convection mechanism. After a review of the 
information available in Tong, it is not entirely clear 
whether the forced convection mechanism or the nucleate 
boiling mechanism is dominant for the relatively low flow 
velocities and low heat flux of this model. In any case, 
the Schrock and Grossman correlation applied at low heat 
flux represents a type of compromise between the entirely 
forced convection assumptions of Dengler and Addoms and 
the "mixed" assumption of Schrock and Grossman. When the 
Schrock and Grossman correlation is applied under the 
conditions of very low heat flux and quality (x < 0.05), 
aS in the first pass of the boiling section of this model, 
1t can be seen that the ratio ee becomes less than 
See This is both unrealistic and computationally diffi- 
cult to handle. Therefore, in the boiling section, the 


model assumes that Dips = h. for the first pass or for any 


Zz 
pass where x < 0.05. This assumption will produce a more 
realistic but still conservative design. 

One other simplification should be explained prior to 


discussion of the actual model calculations for the boiling 


section. Since the boiling section calculations are performed 
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for one complete pass at a time, the average quality, Kos oat 
a particular pass 1s used as the "local" quality to compute 
the inside heat transfer coefficient from the Schrock and 
Grossman correlation. 

The first pass in the boiling section, as indicated 
in the calculations for the heating section, will, in 
general, involve both heating water to the saturation tem- 
perature and the initiation of boiling. From the calcula- 
tions for the last pass of the heating section we have the 
following quantities for the first pass of the boiling 
section: water inlet temperature We. 
outlet temperature (rt ). The gas side heat transfer 


out 
coefficient, Bet 1s obtained, using the average gas tem- 


), average gas 


perature in the boiling section, in the same way as for 
the heating section. 

A rough estimate for the gas temperature into the first 
pass must be obtained in order to begin calculations for 
that pass. This is accomplished by calculating a rough 
overall heat transfer coefficient, a4 for the pass using 
the tube wall resistance, Les and neglecting inside resis- 
tance. A pass effectiveness, co is then calculated assuming 
that only boiling takes place in the pass. The expression 
for effectiveness for fluids unmixed between passes for 


C = 0 becomes 


E (1 - exp (-N) 
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where N = NTU. Using ED calculated from this expression, 


the rough gas temperature for the first pass is 


eae a 


rt «op fin ~ "Sout 
Jin ee 
p 


Since this pass has been simplified by allowing the inside 
heat transfer coefficient to be computed by the Dittus- 


Boelter correlation this initial guess for tt i Seeai@ite 


essential. It would, however, become Area ane 1f we were 
attempting to "Split" the pass with two inside heat transfer 
coefficients, one for heating and one for boiling. InStead, 
the inside heat transfer coefficient, er is calculated for 
the entire pass using the Dittus-Boelter correlation with 


water conditions corresponding to the average water condi- 


tions in the heating section of the pass, that is, 


With this inside heat transfer coefficient, an expression 
fer the overall heat transfer coefficient for the first 


pass may be written, 





yt ~ il 
Ol l . A. ry a ds . Bip 1 
ne Dy ie Ne 1D on ni 


The pass effectiveness, ae is now calculated in the same 


Manner aS previously described. With ene the average gas 


45 





temperature in 1s calculated from 


ML iL 3h 
pio =p Smin “fin 7 “max “Sout 
Jin et i = 


e 
Pp “min max 


The pass heat transfer may be calculated from 


Pmewouclet quality is calculated from 


zt 2. a . ii hayek 
out h : 
fg 


In order to determine the location of the interface between 
heating and boiling the total inside area devoted to heating 


in the first pass is calculated from 


it 
i Neatt ~ fin 
A e oon. ( ) A e 
aim 1 1p 
h h 
se se 
Suc 1B et 


Subsequent fluid passes in the boiling section are 
calculated until the quality at the outlet of a pass is 
found to be greater than 1.0. When this occurs, the boiling 
section is allowed to end with the last pass for which 


Xout < 1-0- A "mixed" pass where both boiling and superheating 
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takes places will be calculated in the superheating section 
of the model. 

After the first pass, boiling section fluid passes are 
calculated in the following manner. 
1. The steam quality at the outlet of the previous pass 
becomes the inlet quality for the current pass, Xen: The 
average gas temperature at the inlet to the previous pass 
becomes the outlet average gas temperature, T pot Oimaeene 

Jout 

SuLrent pass. 
2. Using the last overall heat transfer coefficient calcu- 
lated for the previous pass, Usa approximate values of 
NTU and ED are calculated. 


3. With this approximate Es the average gas temperature 


into the current pass is calculated from 


€ TT. ae 
T = P Sac Jout 
Jin =D ee 
4. Using this gas temperature in, T , the pass heat 


niet 
transfer, Qo is calculated. This allows the calculation 


Prana aval! cy Out, x , aS well as 


of enthalpy out, h Gut 


out 


the average quality in the current pass, x 


5. The inside heat transfer coefficient, h may now be 


tpt’ 


calculated, 


2 (as 1g 
hoe = Bla + Als) bh, , 
P fg tt 
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where the Martinelli parameter oot 46 Calculated using 
average quality, xo etigwecel. dade terms are as previously 
described. 
6. A new overall heat transfer coefficient is now calcu- 
lated from which a new pass effectiveness can be calculated. 
7. This new pass effectiveness yields a new average gas 
temperature into the pass from which a new pass heat 
transfer may be calculated. 
8. The new pass heat transfer is compared with the pre- 
viously calculated pass heat transfer. If the two heat 
transfer calculations do not agree to within 5% the entire 
set of calculations is repeated. 
9. When the pass heat transfer calculation converges, the 
final quality at the outlet of the pass is calculated. If 
the quality is less than 1.0 another pass is calculated for 
the boiling section. 
10. If the quality is 1.0, the boiling section ends with 
the current pass. If the quality out exceeds 1.0 the boiling 
section ends with the previous pass. 

Calculations for the average gas film temperature and 
gas-side pressure drop are performed in the same manner as 


for the heating section. 


mee SUPER RHEATING SECTION 
The first pass of the superheating section will, in 
general, involve both boiling and superheating. The inside 


heat transfer coefficients will be quite different for the 
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fouling and Superheating portions for the pass, and, for 

this reason, the pass will be "split" to allow more accurate 

calculation of the heat transfer in these two regions. An 

Outline of the calculation scheme is as follows: 

1. Calculate the gas-side heat transfer coefficient for 

the section in the same manner as for the previous two 

sections. 

2. Estimate the average gas temperature into the first 

pass. This is only a rough guess for T 4 to begin the 
in 

ea culation. 


3. Calculate the required heat transfer in the boiling 


portion 


where he 1s the enthalpy out of the last pass of the 
a 


boiling section and he is the enthalpy of saturated 


atv 
mia DOL . 

4. Using ome and the average quality in the boiling portion 
of the first pass, calculate the inside heat transfer 


Seetficient, h uSing the Schrock and Grossman correla- 


wena’ 
tion, aS in the boiling section of the model. This, together 
with the gas-side heat transfer coefficient and the wall 
resistance, will yield an overall heat transfer coefficient 
Beuetne bolling section, oe 

5. The calculation of the area required for the boiling 


Portion of the pass is performed in an iterative scheme as follows. 
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a. An initial estimate for the required boiling area, 
Al, is made. 
b. The mass flow rate of the gas over the boiling 


portion is then 


c. An effectiveness for the boiling portion may be 


calculated from 


: = rb 
b b 
ear - T ) 
Be Jin wea 
where 
2 2 & ee 
min Gabg 


d. Using the previously calculated overall heat 
transfer coefficient, the estimated boiling area and the 


consequent gas-Side heat capacity, NTU is calculated from 


e. A "test" effectiveness, en can be found by applying 


this NTU in the relationship for effectiveness 


Ey = eee Or! On. 


: : ic 
f. The two effectiveness calculations, Eb and Ens are 


compared. 
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(10) Ae te oe 1s greater than Eps 


estimate for the boiling area, Al, was too large for the 


then the original 


heat transfer required and the specified gas temperature 

into the pass. So, the estimate for the boiling area is 
decreased and a new Al is established. 

= 

(2) Lt EL 5! 


g. A new gas-side heat capacity is calculated from 


1s less than ¢€ Al is increased. 


Cae a ih 
h. A new boiling portion effectiveness is calculated 


uSing the latest boiling area estimate 


soe sat 
1. As before, an NTU is calculated with the new Al. 
This generates a new "test" effectiveness, a, which is 


; ic 
compared with ¢€ Tf ¢€, and €, are not equal, a new Al is 


is) b b 
calculated based on the comparison. 
j- These calculations are continued until orn = Ep: 
The current Al is then set equal to ae 
6. Once the area required for boiling, Ae, Ponsa sod teu bar 
average gas inlet temperature has been established, the 


following first pass quantities can be calculated: 


a. A. = area available for superheating 


Syl 





b. m = gas flow rate over the boiling portion 


gas flow rate over the superheating portion 


o 


ale Tt? =) gas temperature out of the boiling portion 
Jout 
me eT iy Zeb 
Jout oe ih, ke 
G bg 
e. yee = gas temperature out of the superheating 
Jout 
POrrelon 
A. St = ae T° 
psh P  Jout Jout 
=) aera San 
Jout Ae 
al 
sh b 
(The weighted average of T lige, aw must equal 
Oulte out 
the T previously calculated.) 


out 
7. Since the heat transfer to the cold side must satisfy 


an energy balance involving the average temperatures of the 
gas, the expected heat transfer in the superheating portion 


may be calculated from the heat release on the gas side 


Ol SS Se eae ORE - T 1 
Pp Ss) ieGh Sha eat 


and the required heat transfer in the superheating portion 


as 
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8. The enthalpy and temperature of the steam out of the 
first pass may be calculated from the assumed heat transfer 
in the superheating portion, Onn: 

9. The steam properties for the superheating portion may 
now be found. These properties, along with bs and wall 
resistance, will yield an overall heat transfer coefficient, 
ae for the superheating portion. The Dittus-Boelter 
correlation is applied to calculate the inside heat trans- 
fer coefficient. 

10. The heat capacities on the gas and fluid sides are 


calculated, using the current gas mass flow rate over the 


superheating portion, fig to calculate the gas heat capacity. 


The heat capacities along with rae and asm yield the effec- 
tiveness for the superheating portion of the pass, E oh’ 

ll. A new average gas temperature into the pass may now 

be calculated from the superheating portion effectiveness 
and the superheating portion heat transfer required to 


Satisfy the original heat balance from the expression 


Q 
* 
T = Peper c) 1s 
Jin ae € csh 
alam olel 


* 
This temperature, T , 1S a measure of the temperature 

in 
which would be required to produce the originally specified 
Qon under the actual conditions of the heat transfer char- 


acteristics of the superheating portion. 
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* 
ie tit the T 1s the same as the originally specified 
aig 
average gas temperature into the pass, T , then the heat 
a) 
Si 


under the current distribution of area between boiling and 


transfer in the superheating portion, Q 1S possible 
Superheating and the heat balance is satisfied between steam 
and gas sides for the specified gas temperature in. [In 

this event, the calculations for the first pass of the 
sSuperheating section are complete. 


* 
ime If T is less/greater than the current T the 
Jin ae) 


average gas temperature into the pass is increased/decreased 


and the entire set of calculations is performed again until 


14. The second and subsequent superheating section passes 
are calculated in much the same way as in the heating section. 
That is, an overall pass effectiveness is calculated for a 
Specified number of passes, uSing the multipass effectiveness 
relationship from Kays and London [Ref. 9]. 

15. The average steam-side heat transfer coefficient, Nes 

is calculated using the Dittus-Boelter correlation with 

steam properties at steam bulk temperature for the remainder 
of the superheating section. The overall heat transfer 
coefficient and the average pass effectiveness, ace are then 
ealculated. 

16. The formulation for overall effectiveness with the 


number of passes, n, variable, is 
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min’ max,n 
ae ) etl 
E = P 
Oa Peer ca, §/C Ce ° 
( Pp. min max, n _ | Gan 
le-e C 
Pp max 


As in the heating section this expression is used in an 
iterative fashion, beginning with n = l. 
#7. Each time xe is solved, the gas temperature into the 


superheating section is calculated from 


oa min Te ~ “max ss 
ap = in Sout 
3 Co 
in oa min max 
where T 1s the temperature of the steam into the second 
T39) 
pass of the superheater and T is the average gas tem- 
ete 


perature out of the second pass of the superheater. This 
gas temperature in is compared with the originally specified 


gas temperature into the heat exchanger, T 


gl” 

18. If the calculated gas temperature in, T or 1s less 
than silk the number of passes 1s increased semen and the 
calculation is performed again. If T | exceeds Hake 
superheater is allowed to stop with er ase pass. 


the 


Calculations for the average gas film temperature and 
the gas-side pressure drop are performed in the same manner 


as for the heating and boiling sections. 


F. MATCHING SUPERHEATER GAS INLET TEMPERATURE 
The model calculations thus far lead to a waste heat 


recovery unit design for which the gas temperature at the 
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inlet is, in general, lower than that specified in the initial 
conditions. This condition occurs because the superheater 
ends with the last complete pass calculated for which the 
entering gas temperature is less than or equal to the 
originally specified WHRU inlet temperature. This outcome 
is, of course, unsatisfactory since the inlet gas temperature 
is an independent quantity which is a function of the horse- 
power setting on the gas turbine. An iterative scheme is 
employed for matching the T | calculated in the superheating 
section with the eal See ead fi the Inttr1al conditions. 
Since the model does not allow the calculation of frac- 
tional passes and since the scheme for the calculation of 
superheater passes will not allow inclusion of a pass which 
would require a gas inlet temperature higher than that 
specified in the initial conditions, the final gas inlet 
temperature calculated for the WHRU will generally be lower 
than that initially specified. In order to match this 
Superheater gas inlet temperature with that specified 
initially, the heat exchanger gas outlet temperature is 
allowed to rise. Each time the gas outlet temperature is 
increased, the entire waste heat recovery unit calculation 
1s performed again, and the calculated superheater gas 
inlet temperature is compared with the gas turbine exhaust 


temperature. When a match is achieved the design is 


fixed. 
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G. OFF-DESIGN CALCULATIONS 

Once a particular design has been selected the designer 
must test that design at several off-design points. A 
design will define the following WHRU characteristics: 

Dimensions (length, width, height) 
Operating pressure 

Superheater outlet temperature 
Piet Moe CONnf guration. 

The off-design calculation requires that the selected 
WHRU design performance be predicted for a gas inlet tempera- 
ture and flow rate different from that which was used to 
produce the selected design. The design model can be used 
in an iterative fashion to make this calculation. All 
model calculations are made as prescribed in the previous 
sections, with the exception of the pinch point calculations 
which are not performed. 

The off-design procedure requires that the designer 
fix the WHRU frontal dimensions, operating pressure, fin- 
tube configuration and minimum gas outlet temperature. A 
design is then produced for a gas inlet temperature and 
Flow rate at an off-design point. This design is checked 
to determine whether it conforms to the height dimension 
of the WHRU designed at the gas conditions of the off-design 
point. If the height dimension is not matched then the 
number of WHRU passes must be adjusted until a match with 


the design-point WHRU is achieved. The number of WHRU passes 
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may be altered by either adjusting the WHRU gas outlet 
temperature or the superheater outlet steam temperature. 

The effect of either of these adjustments is to change 

the steam flow rate and the heat transfer rate. It would 

be possible to adjust the steam flow rate and allow the 
superheater steam outlet temperature to remain fixed. 
However, for the method of this model, the computation time 
required for this means of control would be large. The 
actual procedure used produces designs relatively quickly 
which closely approximate the results which would be 
Obtained by a fine control on steam flow. Each time the 
Superheater steam outlet or gas outlet temperature is adjusted, 
a new design is produced. This procedure is continued 

until a design is produced for which the height (number of 
passes) equals that of the deSign point WHRU. When the WHRU 
height at the off-design gas conditions is equal to that 

of the design point WHRU, the gas inlet temperature is 
matched with the gas inlet temperature specified for the 


BEE-design point, using the procedure described in Section E. 


ree COGAS SYSTEM OUTPUT MODEL 

The COGAS system output model calculates performance 
parameters for the combined steam and gas turbine system for 
a specified gas turbine input power. The WHRU model pro- 
vides a steam flow rate, pressure, and temperature for a 
Specified set of gas turbine exhaust gas conditions which 


correspond to a particular gas turbine horsepower setting. 
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This WHRU output is applied in a simple Rankine cycle to 
calculate the steam-side power output. The gas-side 
Beessure drop calculated in the WHRU model is used to 
arrive at a revised gas turbine horsepower. By combining 
these two power outputs, a COGAS system power can be calcu- 
Meeea. Other performance indicators such as specific fuel 
consumption, thermal efficiency and steam turbine share of 
the load are also calculated in the COGAS system output 


model. 


1. Rankine Cycle 


The following conditions are assumed for the Rankine 
cycle calculations: 


a. No Line losses 


b. Steam turbine efficiency = es C535 
e. pendenser pressure = 2.0 psia (4 in Hg) 
d. Feedwater heater pressure = 15 psia 

e. Feedwater temperature = 200 F 


f. Pumping power required is neglected 
Gg. Fuel LHV = 18400 Btu/lbm 

ihewtolilewing input conditions are received from the 
WHRU model: steam turbine inlet pressure (Pi), steam turbine 
inlet temperature (Ti), steam flow rate (Me), gas turbine 


exhaust pressure drop Ve at From Py and Ty the turbine 


inlet enthalpy (h and inlet entropy (s,) can be found. 


1? i 
The turbine exhaust steam quality, assuming isentropic 


expansion, x is calculated from 


vic” 
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where 


i 


S¢ entropy of saturated water 


Sig entropy increment for evaporation. 


The steam turbine exhaust enthalpy assuming isentropic 


expansion, h may now be calculated from 


2s’ 
Nos ~ Ne 3 Xs Neg 
where 
ne = enthalpy of saturated water 
Be = enthalpy increment for evaporation. 


The isentropic turbine work is 
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and 


In this model it is assumed that heating steam for the 
feed-water heater is provided, via a reducing station, from 
the main steam line. The fraction of the steam mass flow 


rate from the WHRU required for the feedwater heater is 


mee fw 2 
h, ~ hy 


where 


2% 
ll 


Fw enthalpy of the feedwater 


DP 
ll 


enthalpy of the condensate at feedwater 
heater pressure. 


Therefore, the steam turbine power output is 


2. Gas Turbine Performance Calculations 
The original gas turbine horsepower, exhaust gas 
flow rate, and exhaust gas temperature inputs to the WHRU 
model were hand-calculated using figure 4 of Ref. 11 and 
Proures 1, 6, and 8 of Ref. 10. The assumptions used to 


calculate the WHRU model inputs were: 
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Ambient temperature = 100 F 
Ambient pressure = 14.696 psia 


Humidity = 0.0 


II 


Fuel LHV 18400 Btu/lbm 


Inlet loss = 4.0 in. HO 


Exhaust loss = 6.0 in. HO 
In the COGAS system output model, the gas turbine horse- 
power correction factor for the gas-side pressure drop in 


the WHRU (BHP for 6" loss/BHP for WHRU loss) is found by 


solving 


Cphhp See O25 + OF002125 Ne ae F 


where chhp = BHP for 6" loss/BHP for WHRU loss. This rela- 


tionship was found by extending the duct loss relationship 
Memeig., 7, Ref. ll linearly from the point of 6 in. HO 
exhaust duct loss already assumed. The final gas turbine 


horsepower is 





* 
where © ai is the input gas turbine horsepower. The total 


COGAS system horsepower, P is found by adding the steam 


cour 
turbine and final gas turbine horsepowers. 
In order to calculate the gas turbine specific fuel 


consumption fig. 4 of Ref. ll is entered with the final 


Gr 





Gas turbine horsepower, and the s.f.c. is read from the 
Peopeller law curve superimposed on that figure. The s.f.c. 
is 


correction factor for exhaust loss in the WHRU, Clear 


found in the same manner as Osis PCOM lg. so Ob she ela), 


where 


Cif See (Sens.c. fOr total exhaust AP)/(s.f.c. for 6" exhaust AP) 


and the gas turbine specific fuel consumption is 


The specific fuel consumption for the combined system is 
calculated by first finding the fuel consumption rate of 


the gas turbine 


Mevel Sol Go ae 5 sae 


The COGAS system s.f.c. is then calculated from 


_ “fuel 
COGAS Paes 





To allow the direct comparison of the COGAS system with the 


all gas turbine system the specific fuel consumption of the 


* 


gas turbine at the new, higher COGAS system power, Sato fet 


is calculated in the same manner as SK ofa ai: 


6 3 





The thermal efficiencies for the gas turbine at 
input power, the COGAS system, and the gas turbine at COGAS 


System power are calculated from the relationship 


ss 2525 
Teh Sis ce itv 


Additionally, the steam turbine share of the load 


is calculated from 2 a 
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Tift. RESULTS AND CONCLUSIONS 


A. BACKGROUND 

The foregoing model was applied in a computer simulation 
program written in FORTRAN IV for the IBM 360-67 computer. 

A listing of this program appears in Appendix A. Addition- 
ally a set of supporting programs was written for water, 
steam, and air properties. A listing of this set of pro- 
grams appears in Appendix B. 

In order to test the model and develop an understanding 
of the behavior of the WHRU, a set of designs were produced 
for a variety of initial conditions. These designs were 
then reviewed, and an analysis of the model behavior was 
performed. After completion of the model analysis, several 
of the designs were selected for observation at off-design 
conditions. The off-design simulations were produced for 
gas turbine input powers corresponding to ship speeds of 
mo, and 20 knots, which represents the entire range of 
gas turbine input powers considered in this thesis. Finaliy, 
two designs were selected for detailed testing at off-design 
conditions. These latter two designs were tested at gas 
turbine input horsepowers corresponding to speeds of 9 to 


20 knots, at one knot increments. 


B. DESIGN VARIABLES 
The following set of design variables was selected to 


produce the basic set of designs. 
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1. WHRU operating pressure: Three WHRU operating pressures 
of 400, 600, and 800 psia were selected for design production. 
Bee WHRU geometric scale: In order to test the model 

response to a range of fin-tube sizes, designs were produced 
for the basic fin-tube configuration described in the geometry 
subsection at 1.0, 0.75, and 0.50 scale. 

3. Frontal dimensions: The approximate WHRU space constraints 
for a DD-963-type engineroom are shown in figures 7, 8 and 

9. Consistent with these constraints, frontal dimensions 
mieel2 x 12° and 12" x 15' were selected for the initial 
design set. 

4. Gas turbine input power: The gas turbine input horse- 
power determines the input gas conditions, temperature and 
Flow rate, for the WHRU. Since the most advantageous use 

of the COGAS system was assumed to be for "cruise" speed 
conditions, the range of gas turbine input horsepower con- 
Sidered corresponds to the gas turbine power required for 

a ship speed range of 9 to 20 knots, with the gas turbine 
Operating alone. The three gas turbine input powers selected 


for the production of WHRU designs are shown in Table I. 


Gas Turbine Input Powers 


GTHP Approx. Speed (kts) Temperature (F) Flow Rate (lbm/hr) 
16421 Ze 849 407589 
8526 16 yAg 328641 
1684 9 689 159731 
Table I 
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SeeeewHRU gas Outlet temperature: The WHRU gas outlet 
temperature is fixed at 400 F. Reference 3 states that 

the principal cause of "hot" corrosion in the WHRU would 

be due to the formation of H,SO, from an expected concentra- 
tion of SO, of up to 50 parts per million. Reference 13 
states that dew point condensation of vapor in the exhaust 
Bes becomes a problem for corrosion if the stack temperature 
is depressed below 250 F. Therefore, a minimum average gas 
outlet temperature of 400 F is considered to be reasonably 
safe from the aspect of gas-Side corrosion prevention. It 
should be further noted that the outlet gas temperature 
calculated by the program iS an average temperature, based 
on the heat release from the gas. Therefore, the actual 
temperature of the gas in the stack will be cooler than 

400 F in some locations due to the non-uniform distribution 
of the gas temperature as it leaves the WHRU. It should 
further be noted that, although this temperature is a fixed 
design variable for design production, it will vary somewhat 
from 400 F in response to the program calculations for 
Peeing the minimum pinch point and for matching the gas 
inlet temperature with that specified in the initial 
@ondiltions. 

6. Superheater outlet steam temperature: The target WHRU 
Superheater outlet steam temperature is 650 F. This tempera- 
ture will also vary somewhat in response to steam flow rate 


changes brought about by the program calculations to match 
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the gas inlet temperature with that specified in the 
iartial conditions. 
7. WHRU water inlet temperature: The WHRU water inlet 
temperature is fixed at 200 F. 
8. The other major COGAS system parameters are fixed as 
follows: 

Condenser pressure: 4 in. Hg 

Steam turbine efficiency: 0.85 

Feedwater heater pressure: 15.0 psia 


Fuel LHV: 18400 Btu/lbm. 


fee tat DESIGN SET 

A set of 54 designs was produced for the combination of 
design variables just discussed. The array of design 
variable combinations is described in Table II. The summary 
computer output pages for the results of each design run 
are presented in Appendix C. 

As discussed in Section A, the purposes for producing 
this design set were to test the model and to promote an 
understanding of the likely behavior of a WHRU for a range 
of conditions. To this end, the following set of output 
variables was considered. 
fm WHRU output 

a. height (number of passes) 

b. gas inlet temperature 

Ge gas outlet temperature 


dd. steam flow rate 


cae 
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e. gas side pressure drop 
f. pinch point temperature difference 
g- heat transfer rate 
2. COGAS system output 
a. gas turbine horsepower (including loss from the 
WHRU gas-side pressure drop) 
b. steam turbine horsepower 
c. COGAS system horsepower 
d. steam turbine share of the load 
e. specific fuel consumption 
(1) s.£.c. of COGAS system 
(2) s.f£.c. of gas turbine at COGAS system horsepower 
ff. thermal efficiency 
CL) ea) of COGAS system 
(2 Neh of gas turbine at COGAS system horsepower 
Clearly, several of these variables are related, such as 
mer. c. and Neh’ paduener additional information 1s previded by 
Studying all variables ina closely related group. Therefore, 
the following subset of output variables was selected for 
analysis. 
Ll. WHRU height 
2. Steam flow rate 
3. Gas-side pressure drop 
fer inch point temperature difference 
9. WHRU heat transfer rate 


6. WHRU gas outlet temperature 
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7. Steam turbine horsepower 
8. COGAS system horsepower 
9. COGAS specific fuel consumption 

The above output variables are presented in tabular form 
using the format of Table II. Table II also serves as an 
index for the set of basic designs. Additionally, the 
Significant trends in the variables are presented graph- 
ically. Where the trend of the results was ascertained to 
be the same for the two frontal dimension sets, only the 
values for the 12' x 12' front were presented graphically. 
Where clarity was better served, only one scale was plotted 
as representative of the trend described. 
meeect of Scaling. 

The only significant trends noted with respect to scaling 
were a decrease in WHRU height with a decrease in scale 
and the same trend for gas side pressure drop (see Tables 
TII and IV). The trends are shown in figures 10 and ll. 

The basic fin-tube dimensions of the model were scaled 
meetactors of 0.75 and 0.50. For any scale with constant 
frontal dimensions, the total inside and outside areas 
remain the same. That is, for a constant length tube, the 
scaled inside/outside area (A? /A2) is proportional to the 
full scale areas (A, /A,) with the scale factor, s, as the 


constant of proportionality. Hence 
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Since constant frontal dimensions are to be maintained, 
the number of tubes per row is scaled up by the scale 


factor, so 





N N 
S . _t/p ae 
iip A S A; and oye = z S A, ; 


where Nie /p 1s the number of tubes per pass. Therefore, the 
scaled outside and inside areas for a pass for constant 
frontal dimensions are the same as the full scale areas. 
Now, the scaled cross-sectional area for fluid flow on the 


inside will be, 


which is equal to the Reynolds number for the unscaled 
geometry. It can also be shown that the outSide Reynolds 
number remains constant with scaling. Recalling that the 
Dittus-Boelter correlation was used to calculate the inside 
heat transfer coefficient for the single phase regions and 
used with a multiplier for the two phase heat transfer 
coefficient and assuming constant fluid properties, the 


Nusselt number remains constant 
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So, the scaled inside heat transfer coefficient 





= aes eee 
wae cede ” s 
i 
increases with a decreasing scale factor. Similarly, for 


the scaled gas-side heat transfer coefficient, 


an increase is seen with a decreasing scale factor. 
Therefore, in addition to a decrease in height attributable 
to a scaling down of the fin-tube dimensions, there is also 
an improvement in heat transfer which results in a further 
reduction in height by elimination of passes. This trend 
is supported by the results. 

Recalling the relationship used to calculate the gas- 


Side pressure drop, 


it can be shown that the scaled gas-side pressure drop is 


inversely proportional to the scale factor, 
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where C iS a constant based on the assumption of constant 
gas properties between scales. For the same number of 
passes, then, an increase in gas-side pressure drop could 
be expected. However, since the number of passes is decreased 
with decreasing scale, a net decrease in gas-side pressure 
drop was experienced (Fig. 11). 

Some differences in pinch point with scale were observed 
in Table V, but no consistent trend could be established. 
It was expected that pinch point AT would remain fairly 
constant with scale for fixed pressure and input gas condi- 
tions. That result was apparent to some extent, particularly 
in the range of high gas temperature and low steam pressure. 
This finding is consistent with the manner in which the pinch 
point AT is established. That is, the program fixes an 
fateial pinch point AT only when the rough pinch point AT 
is less than 25 F in the initial calculations. At the lower 
pressures and higher gas temperatures this mechanism would 
not be operative. Also, the final and initial pinch point 
AT calculations are based on the average gas temperatures 
at the inlet to the pass where the AT is calculated. There- 
Fore, the actual pinch point may be either lower or higher 
than the one calculated, depending on the actual gas temp- 
erature distribution at the inlet to the pass. The pinch 
point AT will be discussed at greater length in the analy- 


Sis of the effects of pressure. 
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Effect of Pressure. 

The most important effect of pressure on the model 
performance arises from the calculations which establish 
the minimum pinch point temperature difference. As related 
in the model description, the minimum pinch point AT is set 
at 25 F. The WHRU outlet gas temperature is adjusted in 
the initial model calculations in order to ensure this 
minimum AT. The initial model calculations establish the 
steam flow rate from a simple energy balance involving the 


gas inlet and outlet temperatures Booed ), the gas flow 


3 


rate, the water inlet temperature (T,,), the required steam 


ial 


outlet temperature (Tey), and the steam pressure. Tentative 


interim gas temperatures for the heater inlet, T and 


Sif 


boiling section inlet, T are also established in this 


G2o 


calculation set. In order to establish the minimum pinch 
point AT, the model tests the interim heater gas inlet 


temperature against the criterion 


where Te is the temperature of saturated water at the 


pressure specified. If nee is less than Te5 + 25 then seg 


is set equal to T + 25 and a new gas outlet temperature 


ie 


1s calculated as follows, 


a5 





where 


an el 
SS 2 hae 
£4 a 


As an aid to understanding the behavior of the model, it 

is useful to predict the WHRU gas inlet temperature at 

which this pinch point calculation becomes operative for 

a particular pressure. This is accomplished by establishing 
a "Critical" heater gas inlet temperature, to for each 


pressure considered, 


where T is the temperature of saturated water at the 


f sat 
pressure specified. The ratio a can also be solved for 
each pressure considered, based on a constant water inlet 


temperature of 200 F and a constant superheater outlet 


temperature of 650 F, 


“a Ae = Ay 
eo on 
4 il 
where hz = enthalpy of saturated water at the pressure under 


i 


consideration. Finally, assuming a minimum WHRU gas outlet 
temperature of 4-0 F the minimum WHRU gas inlet temperature 
at which the pinch point calculation comes into effect may 


be calculated as follows, 
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At - (1 - a) 400 
gl O 


Table VI summarizes these calculations for each of the pres- 


sures considered in the design set. 


Minimum Gas Inlet Tenweeramure fOr Pinch Point Calewlation 


* * 


Pressure ft GE) a Al Cs) 
~(psia) a % : 
400 470 5 ZAAGS, 720 
600 497 7262 Ae 
800 543 7209 878 
Table VI 


Now the pressure at which the pinch point calculation becomes 
operative for each gas turbine input power considered is 
predicted. The pressure for which the pinch point calcula- 
tion begins to adjust gas outlet temperature is listed below 


for each gas turbine input power considered, in Table VII. 


Pressure At Which Pinch Point Calculation Takes Effect 


Gas Turbine AE (F) Pressure T (EF) 
Input Power gt (psia) SD 
16421 849 800 878 
8526 742 600 770 
1684 689 400 720 

Table VII 
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It should be recalled that the outlet gas temperature is 

also adjusted by the iterative technique employed to match 
the model-calculated gas inlet temperature with the gas 

inlet temperature of the initial conditions, once the 

design is established. This is a relatively minor adjust- 
ment, however, and 1s never greater than 10 F. These adjust- 
ments to WHRU gas outlet temperature are apparent in the 
values of Table VIII and are displayed graphically for the 
/meeex 12° frontal dimentions and 0.75 scale in Fig. 12. 

The final calculated pinch point for each of the design 
sets reflected the influence of the initial pinch point 
calculation. These results were shown in Table V and Figure 
13. It is re-emphasized here that the calculated pinch 
point AT is useful only as an indicator of what the actual 
pinch point might be for a particular heat exchanger design. 
This 1S true because the model neglects the possibility of 
a gas temperature distribution across the inlet of a WHRU 
pass and, instead, uses the average gas temperature for 
ea culations. 

There are two other trends in the design results which 
emanated directly from the control exercised by the model 
Over the pinch point AT and the resulting gas-fluid tempera- 
ture distributions in the deSigns produced. These trends 
were apparent when the effects of pressure on WHRU height 


and WHRU heat transfer/steam flow rate were examined. 
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The effects of pressure on WHRU heat transfer and steam 
flow rate are shown in Figures 14 and 15 respectively. 
Tables IX and X list the actual values obtained for WHRU 
heat transfer and steam flow rate for the designs produced. 
WHRU heat transfer and steam flow rate directly reflected 
the adjustment of the WHRU gas outlet temperature to main- 
tain the minimum pinch point AT. That is, as the WHRU 
outlet gas temperature was raised, the heat transfer rate 
and steam flow decreased. The small increase shown in 
Figures 14 and 15 in heat transfer and steam flow rate 
for 16421 gas turbine input power at the 0.50 scale reflects 
only a minor adjustment difference in T for matching 
v- between 400 psia and 600 psia. a 

the relationship of WHRU height to WHRU pressure is 
Slightly more complex. Figures 16, 17, and 18 demonstrate 
the approximate gas-fluid temperature distribution along 
the total length of the WHRU for the 0.75 scale at the 
three pressures and input horsepowers considered. The 
general trend observed in these diagrams of decreasing gas- 
Fluid temperature difference with increasing pressure is 
representative of all the designs produced. This decrease 
in gas-fluid temperature difference, taken by itself, would 
produce larger heat exchangers with increasing pressure. 
This trend may be observed for the 12' x 12' front at the 


high gas turbine input horsepower in Figure 19. At the 


medium and low gas turbine horsepowers, however, the WHRU 
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height remained fairly constant or actually decreased with 
increaSing pressure. This trend is explained by the 
decreasing heat transfer rate with increasing pressure 
attributable to the pinch point calculation adjustment to 
WHRU gas outlet temperature shown in Figure 12. 

The design set results for steam turbine horsepower 
output are given in Table XI. These results are presented 
@eadpnically for the 12° x 12° front in Figure 20. Two 
separate trends were identified. The steam turbine horse- 
power demonstrated an overall increase with increasing 
pressure for the high gas turbine input horsepower. Although 
there is a small net decrease in WHRU heat transfer rate 
with increasing pressure at the high gas turbine input 
power, the steam is rejected to the condenser at a lower 
quality for the higher pressures. Therefore, there is less 
heat rejection in the condenser for higher pressures. This 
better performance of the high pressure systems in the 
Rankine cycle more than compensates for the small decrease 
in WHRU output at these pressures. -For the medium and low 
gas turbine input powers the trend of decreasing WHRU output 
with increasing pressure (Fig. 20) was dominant in the 
Rankine cycle also, producing systems of decreasing steam 
turbine horsepower with increaSing pressure. 

The design set results for the COGAS system combined 
horsepower are given in Table XII. The results are displayed 


Metpatcally in Figure 21 for the 12' x 12' front. It is 


99 











876 SST 

S501 CLOT 
yer ec iae 
LOLC Livy 
pI8c GES 
bpt6c SD EG 
Dees S6TS 
bpSTS 89TS 
C88P cCBP 

OG O-91Te9S. GZ O=eT129S 
(ST X «oT + FUOTY 


676 
c90T 


p8tt 


OL9C 
T98C 


LS6¢ 


StS 
660S 


678P 





GI E=2T29S 


YAHMOddSUOH 


[IX HIV 


bS6 
8cOTt 


Bish 


ECIC 


8C8BC 


pS6¢ 


Lacs 
GEGs 


Se Sg ba 


0G O=9Te9S 


ie. L.. BX 


Lew 


C96 
OvOT 


C9TCT 


OG9¢ 


BEBC 


EGC 


OOTS 


[80S 


SO8P 


GL°O=9LTP9S 


ANTGYNGL WVALS 


SHUG 


SS6 
TLOT 


CLUtT 


USC 
CIBC 


Voc 


SO0¢cS 
CLOGS 


Geav 


O° T=aTeOS 


008 


LON: 


OOV 


OO8 


SIONS, 


OOV 


O08 


1013) 


OOV 


2a See 


ies 9) IE 


9¢S8 


Tey 


dHbO 


100 





Oc FAUNA 


(OtSad) Saeissacel 


oO rep) aN rm 
8 8 8 S 
a a BBS 
dH esl K—- — — _ ees BBB) 
€ —— =e 
BBcl 
BBBe 
98Se L 
dH 8250 MK sam © 
ge eens. AJ 
—— Ter BBBE W) 
= 7 
a 
3WOS @S°S@ + BBSE © 
aWwOS GZ-a ow oe 
JWOS @a°t x BBBr AW 
Tw ‘2190S 
,at%,e3 *LNOS BBSPb 
ae 2228S 
dH Tepe == —~ _xX- 
BBSsS 


Ae > Sadie) Sie SeIM@elad S10 


Om 





102 





TIX HWIdVva 
Croc 819d IIL 919C BCI OTIC O08 
(eds. Vere VCLC 069C LOC CELTS O09 psot 
GEC CV8c OV BIC 6£8C CCEC CEBC OO” 
860TT 6cOTl 6P0TT 6660T 9660T ETOTT O08 
SOCTII 8IcIl 6tecIl TOCTT LOCTTI TIctl O09 9¢S8 
9OETT TSt tl 2) 2 a 1B SGeUL 96cIT L6cTT OOV 
8IETC EViGiC SOCTC LOCTC 9TTIT¢ COTTC OO8 
CLOLC 69cIC¢ T8 Td 98cTC CTLT¢ 990TC O09 Lev t 
'l60c 67602 99602 TE80C VL80C T980¢ OOV 
OSes ten Se O-2leoS CO 1=]e129S OS O=2T8 °F SZ OF=21Te 4S ~O =e els esse eae dHLO 
ce ee, Ge eA | - JUOT Y me | hens sO 


YAMOdHSHYOH WHLSAS SWDOO 





ram hope iG 


Cyl Sa) weiaiiaosad 


0) 
GQ 
G 


B08 
Gap 
Wese 


ee I oe. oe ame ie ee. SU a e7 SC Spc ee. ,  e 


SHLD B97 
BB02 


BEOb 

8889 

8088 
FHLS 8268 99091 


BBGcI 


oe 





3 wos a8°o + BRBrI 
SWS &2°SO 

SWS SB°t xX 9en91 
Tw '%'sFwWosS 

2O3%,27 PLNOM4 90881 


ce 


* 
| 
| 
| 
| 
| 
4 
| 
| 
| 
| 
| 
| 
k 
Sah CaesaOn 


dHLD T2693 BOBBe 


gaae¢2 
Geet SNe) MG) cS SO i) 0) 


OG 





Meearent that the system horsepower output trend followed 
that for steam turbine horsepower output. This was the 
expected result. This trend is modified only slightly by 
additional gas turbine horsepower losses due to a changing 
gas-Side pressure drop across the WHRU. 

Finally, the system specific fuel Poneunemiion followed 
closely the results for steam turbine horsepower output. 
That is, Slightly better systems were observed at lower 
pressures for the medium and low gas turbine input powers 
and a slightly better system was produced at higher pressures 
for the high gas turbine input horsepower. These trends 
are shown in Figure 22 and Table XIII. 

Effect of Gas Turbine Input Power. 

An increase in WHRU height with increasing gas turbine 
input power was noted in the results for the design set 
produced (see Fig. 23 and Table III). This increase in 
WHRU height generally followed the increase in WHRU heat 
transfer rate with increasing gas turbine input power, Fig. 
24. This trend is modified only by a slightly increasing 
outside heat transfer coefficient with increasing gas turbine 
input power (see calculations for gas-side heat transfer 
coefficient in the model description). 

As expected, the COGAS system specific fuel consumption 
improved (Fig. 25) with increasing gas turbine input power. 
This improvement represents the trend of improvement in 


the gas turbine itself with increasing power, since the 
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contribution of the steam system decreases with increasing 
power. 
Effect of Frontal Dimensions. 

The only significant trends for changing frontal dimen- 
sions were in WHRU height and volume (see Tables III and 
XIV). The expected result of decreased height with increased 
frontal area was observed. The gas-side pressure drop 
follows this decrease in height. The WHRU volume showed a 
Significant increase for increased frontal area for all 
designs produced. This result is consistent with the model 
formulation. That is, as the frontal area is increased 
and the tube length is held constant the area for fluid flow 
is increased with a consequent decrease in Reynolds number 
and inside heat transfer coefficient. The same result is 
seen in the outside heat transfer coefficient. As the 
frontal area is increased, the minimum flow area for the 
gas 1S increased, causing a drop in the gas side Reynolds 


number and heat transfer coefficient. 


fee OF F-DESIGN RESULTS 

The off-design results were produced using the same 
computer simulation program used to produce the design set. 
The method used for off-design runs was described in Section 
F of the model description. As the computer program is 
presently designed, the designer must manually adjust the 
inputs to produce an off-design run. The procedure is 


reviewed briefly below. 


Alle, 
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First, the designer selects a WHRU design to be inves- 
tigated at off-design conditions. Once this design has 
been selected, the pressure, frontal dimensions, height and 
scale of the WHRU are fixed. The designer then enters the 
gas flow rate and gas inlet temperature corresponding to 
the off-design point to be investigated. A design is pro- 
duced with these gas conditions and the physical charac- 
teristics of the WHRU design under consideration. The 
designer then checks the number of passes (height) for the 
resulting design. If the number of passes does not match 
that of the selected design, the WHRU gas outlet temperature 
and/or superheater steam outlet temperature are adjusted 
to increase or decrease the number of passes until a match 
is achieved. Once the physical characteristics of the off- 
design point WHRU match those of the design-point WHRU, the 
gas inlet temperature is matched with the initial conditions, 
and the performance of the selected WHRU design is established 
for the off-design point under consideration. 

In order to test the feasibility of the procedure described 
above and to demonstrate the off-design performance of a 
wide variety of WHRU designs, ten WHRU designs were selected 
for off-design runs. These designs are listed in Table XVI. 
None of the full-scale designs was considered because the 
heights of the designs produced at that scale were, for the 
most part, either too large or borderline according to the 


Space limitation described in Figure 9. 


iv 





The results of these off-design runs are given in 
Table XVII. The summary output pages for each of these 
off-design runs along with the design-point runs are pro- 
vided in Appendix D. In terms of specific fuel consumption, 
the performance difference between the design sets con- 
sidered was not large. The performance range for fuel 
consumption is even less significant. Table XV demonstrates 


the performance range for the designs considered. 


WHRU Off-Design Performance Ranges 


nicl consump = approx. L£uelaconsump— 
Power Sera. ekeve Cb ejuveou ay, tion difference (gal/hr) 
high 29 =). 340 S792. 6°="7 09.3 Cry, 
medium .391 - .403 4429.2 - 4439.7 ee) 
low SINS) So eopeite: imo4.6 = 1767.5 0.4 

Table XV 


Since the fuel consumption comparison showed very little 
difference between the ten designs considered and since all 
designs considered satisfied the dimensional constraints 
of figures 7, 8, and 9, only two designs were selected for 
additional off-design runs. 

For these additional off-design runs, a pressure of 
600 psia was selected. This pressure did not produce the 
smallest designs but it is a steam pressure with which the 
Navy has some experience in terms of materials requirements, 


Maintenance, and operation. A scale of 0.75 was selected 


LY 





WHRU DESIGNS CONSIDERED 


Design Design 
Run # Pressure =cale Bronte Power 
3) 800 B/S. Ibi LZ med 
14 600 Se: Zac ie med 

16 600 520 ee high 
Z2 400 275 ce hog 
24 400 wi 2 ez sc 1ez low 
SD 800 750 ieZecl> med 
40 600 AN Ls hes lus high 
41 600 ui 2 med 
3)8) 400 0 2s 5 med 

TABLE XVI 
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which yields a tube outside diameter of 1.5 inches and a 

fin Spacing of about 8 fins/inch. The rationale for this 
selection was that a scale of 0.50 with 1.0 inch OD tubes 
and a fin spacing of 12 fins/inch would probably be too 
susceptible to both inside and outside fouling. In order 

to achieve the closest possible comparison, the design-point 
for both designs selected was for the gas conditions corres- 
ponding to the medium (8526 BHP) gas turbine input horse- 
power. Both frontal areas were used for these detailed 
off-design runs. Table XVIII summarizes the characteristics 


of the two designs selected. 


Final Design Characteristics 


design Cie inp 
run # design point scale MeO pressure height 
14 8526 Ney 5 LA SN? 600 Se 
41 8526 Ona75 dls Sram bass 600 29 


Table XVIII 


Table XIX summarizes the performance results for the 
COGAS system with the 12' x 12' front WHRU and Table XxX 
gives the results for the system with the 12' x 15' front 
WHRU. In both cases, the performance characteristics of 
the gas turbine alone, at the COGAS system horsepower, are 
provided for comparison. The summary output pages for each 


of these runs are provided in Appendix E. 
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The specific fuel consumption curves vs. brake horse- 
Mewer for the 12° x 12° front and the 12° x 15' front are 
plotted in figures 26 and 27 respectively. The difference 
in s.f.c. performance between the two WHRU designs con- 
sidered was not Significant. Geometric considerations were 
therefore used in selecting the better of these two deSigns. 
A geometric comparison of the two WHRU designs is given in 


Table XXI. 


Geometric Comparison of the Final Two Designs 


Pic (Orgueelle To@etis Out 
Front Area Height Volume Side Area 
me x15' 181.9 229 528.96 49586.4 
me x12' 144.8 Baz 464.64 42776.5 
Table XXI 


It was concluded that, since both designs met the height 
Menstraint of figure 9, the 12‘ x 12' front was the 

better choice because of the significantly lower volume and 
total outside area requirements. 

After the selection of the 12' x 12' front WHRU as the 
final design, estimates were made of the possible COGAS fuel 
Savings for a DD-963-type destroyer over 1000 hours of 
operation. The NAVSEC Standard Destroyer Profile (Fig. 28) 
was used to determine the number of hours a destroyer-type 


ship spends at each speed during 1000 hours operation. 
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Figure 29 was used to obtain the shaft horsepower require- 
ments for each speed. Finally, figure 30 was used to obtain 
the s.f.c. of the gas turbine operating alone. 

Two estimates of fuel savings were made. The first 
estimate assumed that the COGAS system would be operated for 
speeds of 5 to 20 knots and would be secured outside that 
range. The second estimate employed the COGAS system for 
speeds of 5 to 23 knots. For the latter operating range, 
the criterion for the high COGAS speed was an imposed maxi- 
mum gas turbine horsepower input of 20000 BHP. Both fuel 
Savings estimates are based on the following assumptions 
and simplifications: 

1. The COGAS mode of operation implies the use of one 
engine (COGAS mode) on one shaft with the other shaft 
dragging. 

2. The gas turbine is assumed to operate at idle (1000 BHP) 
at speeds below 8 knots. 

3. Maneuvering combinations, when two main engines would 
normally be on the line, are not considered. 

4. In the speed range of 21 to 27 knots, in the pure gas 
turbine mode, one engine per shaft is used, with two shafts 
on the line. 

5. From 28 to 32 knots four gas turbines are on the line. 

Table XXII provides the estimates of s.f.c. and fuel 
consumption for COGAS between 5 and 20 knots and for pure 


gas turbine for all speeds. The total estimated fuel savings 
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over 1000 hours of operation was 12% which is roughly 
equivalent to 96000 gallons of distillate fuel. 

The estimates of s.f.c. and fuel consumption for pure 
Gas turbine and for COGAS for the 5 to 23 knot speed range 
are given in Table XXIII. The total estimated 1000-hour 
fuel savings for this COGAS operating range was 19% which 
translates to approximately 152000 gallons of distillate 


fuel. 


Bee CONCLUSIONS 

The waste heat recovery unit model and the accompanying 
Simple COGAS system output model provide the basic frame- 
work for additional studies of the application of the COGAS 
system to U. S. Navy gas turbine powered ships. The model, 
as presently formulated, provides a reasonable estimate of 
WHRU size required for the fin-tube configuration considered. 

Estimates of the COGAS system performance in the 5 to 
20 knot range and the 5 to 23 knot range are encouraging. 
Even considering the assumptions involved in making the 
fuel savings estimates, it is likely that a COGAS version 
of a DD-963-type ship would consume on the order of 300000 
to 450000 gallons less fuel in one year of operation than 
Mae current DD-963 class ships. 

The most probable COGAS system would be one which would 
be designed to operate in non-maneuvering situations for 
either speed range considered above. That is, a system which 


would be employed for maximum fuel savings in the cruise mode. 
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During maneuvering and high speed operations the gas 
turbines would be operated alone with the waste heat 
recovery unit dry. Because of the relatively low gas 
turbine exhaust temperatures involved, it is not expected 


that damage would occur in the WHRU when operated dry. 
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IV. RECOMMENDATIONS FOR FURTHER RESEARCH 


A. WHRU DESIGN MODEL IMPROVEMENT AND EXPANSION 

The type of control imposed on the WHRU design model of 
this thesis was essentially the adjustment of the steam flow 
rate to maintain a superheater steam outlet temperature of 
650 F at the pressure specified. This steam flow rate 
control was further extended to accommodate the maintenance 
eea minimum pinch point AT of 25 F. Further investigation 
of WHRU design should include the effects of allowing the 
Superheater outlet steam temperature and/or pressure to 
"float" with changing gas turbine exhaust conditions. The 
combination of controls on the superheater outlet pressure, 
temperature, and flow rate could lead to improved WHRU output 
while still maintaining an acceptable minimum pinch point 
temperature difference. 

In order to increase the precision of the model and lead 
to increased knowledge of the gas temperature distribution 
in the WHRU, the present model could be modified for a "fine 
mesh" approach. That is, instead of solving one pass at a 
time, each pass could be divided into a number of segments. 
The solution could then proceed uSing either the finite 
difference or finite element method. 

Other possible model improvements are as follows: 

1. Include the calculation of fluid/steam-side pressure 


SOD. 
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2. Additional fin-tube configurations should be considered 
for possible improvement of the heat transfer characteristics 
and the weight and space requirements of the WHRU. 

3. Consideration of inside and outside heat transfer sur- 
face fouling should be given when investigating possible 
fin-tube configurations. 

4. Working fluids other than water should be investigated 
for possible enhanced thermodynamic characteristics at a 


minimum cost in system maintainability. 


B. COGAS MODEL IMPROVEMENTS 

After the waste heat recovery unit, the steam turbine 
and condenser should have the highest priority for modeling. 
Both of these units were modeled as "black boxes" in this 
model. The modeling or performance mapping of "state of the 
art" steam turbines for application in the COGAS model would 
allow the designer to further refine the estimate of system 
performance for various WHRU outputs of steam temperature, 
pressure, and flow rate. The inclusion of a steam turbine 
model would also allow the designer to make size and weight 
estimates for that component. The same benefits would be 
derived from the inclusion of a condenser model in the COGAS 
system. As the condenser pressure is reduced, the system 
performance should improve. The improvement must, however, 
be weighed against the increased condenser size necessary 


at the reduced pressure for constant steam and cooling water 


1 ESYS 





conditions. A condenser model would allow the designer 
to make these comparisons. 

In the model of this thesis, pumping power for the 
condensate and feedwater was assumed to be negligible. In 
the case of the condensate pump, this is probably not a bad 
assumption since it is likely that electric pumps could be 
used, and the cost to the ship's electrical output would 
not be large. The feedwater pump, however, would probably 
be steam turbine driven at some cost to the COGAS system 
Output. The required feed pump pumping power for various 
loads should be mapped and included in the model. 

Finally, the likely engineroom machinery layout for the 
COGAS system should be described in enough detail that the 
lengths of piping runs, number of turns, and location of 
valves could be estimated. The inclusion of this feature 
in the model would allow the designer to address pressure 
drop and heat losses between components. This feature would 
also contribute to the COGAS system weight and space predic- 


fon. 


fe SCLSoLEM OPTIMIZATION 

Since the COGAS system is operated with a relatively 
small heat source, optimization should be attempted in the 
design of the system. The technique of non-linear programming 
could be applied to a COGAS system optimization where an 
important system output parameter, say specific fuel con- 


Sumption, 1s expressed as a function of the system design 
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variables such as component weight, component dimensions, 
and gas turbine backpressure. This function is called the 
objective function, and it is maximized or minimized subject 
to constraints which are also expressed as functions of the 
design variables. These constraint functions may be either 
linear or non-linear and express some parameter limit not 

to be exceeded, such as WHRU or condenser total volume. 

For a system of the complexity of the COGAS system a 
"local" optimization approach would probably have to be 
adopted. That is, each of the major sub-models (WHRU, 
steam turbine, and condenser) could be optimized using local 
constraints. Then the entire system could be optimized using 
linking variables, such as enthalpy, flow rates, and tempera- 
tures between components. The process would be repeated 
as many times as necessary to achieve the final system 
optimization. One available computer-based system using 
the technique of non-linear programming is the COPES/CONMIN 


program [{Ref. 14]. 
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APPENDIX A 
MAIN PROGRAM LISTING 
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APPENDIX B 


This appendix contains the program listings for the 
subroutines called by the WHRU deSign program and the COGAS 
model to obtain steam, water, and air properties. Refer- 
ences 15 and 16 provide relationships for use in calcula- 
ting the thermodynamic properties of steam and water. The 
relationships of Ref. 15 for the thermodynamic properties 
of steam were used in Subroutine SS of this appendix. 

The remainder of the subroutines for steam, water and 
alr properties were obtained by fitting polynomials to the 
values for those properties found in References 17 and 18. 
In all cases, the polynomial coefficients were obtained 
using a "canned" program for least squares regression pro- 
vided with the Hewlett-Packard System 9845 Computer installed 
in the Naval Postgraduate School Mechanical Engineering 


Department. 
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